NPS  ARCHIVE 
1969 
READ,  D. 


AN   ITERATIVE   PROCESS   TO  SOLVE 
THE   GRAPH-COLORING   PROBLEM 


by 
Dennis  Spencer  Read 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  93943-5101 


United    Stat 
Naval  Postgraduate  School 


THESIS 


AN  ITERATIVE  PROCESS 

TO 

SOLVE  THE  (3*APH-C0L0RING  PROBLEM 


by 


Dennis  Spencer  Read 


773^£gf 


October  1969 


Tku  document  Ivu  bien  apptuovzd  fri  public  it- 
lexue.  and  &alz;  aM,  distribution  lb  untunitzd. 


o-r>v 


An  Iterative  Process 

to 

Solve  the  Graph-Coloring  Problem 


by 


Dennis  Spencer  JRead 
Lieutenant  Commander,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1959 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
October  1969 


ABSTRACT 

The  intent  of  this  paper  is  to  describe  a  method  of  coloring  a 
map  and  to  present  an  algorithm  for  the  solution  of  this  problem. 
A  computer  program  was  developed  to  provide  solutions  to  the  problem 
of  coloring  a  map  which  consists  of  a  finite  number  of  areas.  This 
algorithm  may  also  be  applied  to  problems  other  than  map-coloring. 
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TABLE  OF  SYMBOLS 


SYMBOL  MEANING 

M  Number  of  equations  used  to  describe  the  border 
constraints  for  a  map. 

N  Number  of  areas  to  be  colored 

x.  Variable  representing  the  color  of  area  i 

A^  The  MxN  matrix  defining  regions  having  common  borders 

a.  The  j—  column  of  AM^ 

a_i  The  i—  row  of  A^ 

A'mxN  ^MxN  +  [aj'  •,,/  aj]lxN 

g .  The  number  of  zero  elements  of  each  column  a .  that 

J  are  not  associated  with  the  constraint  equations 

a .  The  number  of  zero  elements  of  each  column  a . 
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I.   INTRODUCTION 

When  coloring  a  geographical  map,  it  is  customary  to  use  different 
colors  for  adjacent  regions  -  where  adjacent  means  a  common  finite 
boundary,  it  has  been  conjectured,  but  not  proved,  that  "For  any  sub- 
division of  the  plane  into  non-overlapping  regions,  it  is  always  pos- 
sible to  mark  the  regions  with  one  of  the  numbers  0,  1,  2,  3,  in  such 
a  way  that  no  two  adjacent  regions  receive  the  same  number." 

The  statement  that  a  map  can  be  colored  with  four  colors  is 
accredited  to  Moebius,  who  first  proposed  it  in  1840.  This  conjecture 
has  never  been  proven;  however,  a  satisfactory  counter-example  has  not 
been  demonstrated  either.  A  proof  for  five-coloring  a  map  may  be  found 
in  [Courant  and  Rabbins  1941]  and  [Oystein  1967] . 

This  thesis  will  use  the  four-coloring  conjecture  without  proof 
and  present  an  algorithm  that  utilizes  an  iterative  routine  to  determine 
what  color  each  region  should  be.  A  similar  iterative  technique  for 
general  integer  programming  problems  is  demonstrated  in  [Greenberg  April 
1969]  and  [Greenberg  May  1969] ,  which  present  the  proofs.  We  are 
interested  in  showing  whether  the  iterative  technique  can  be  used  for 
large-scale  problems  such  as  the  map-coloring  problem. 

The  problem  formulation  suggested  by  Gomory  [Dantzig  1960]  is  used. 
Let  the  map  regions  be  denoted  by  i  =  1,2, ...  ,N  and  let  x.  be  an  integer 


valued  variable  such  that 


0  i  x,  <  3  (1) 


Courant,  Richard  and  Robbins,  Herbert;  What  is  Mathematics?,  p.  247, 
Oxford  University  Press,  1941. 


where  the  four  values  x.  =  0,  1,  2,  or  3  correspond  to  four  different 
colors.  Since  any  two  adjacent  areas  must  have  different  colors,  the 
constraint  may  be  written  in  an  either/or  form 

either  x.-x.^l        or  x.  -  x.  £.  1         (2) 
if  i  and  j  have  a  common  border.  Equation  (2)  may  be  rewritten  as: 

x±-  x^z  1  -   46±j   (6....  =  0,1) 

(3) 

x.  -  x.  $.-3  +  46. 
3         i         ID 

Instead  of  using  (2)  and  (3)  we  can  simplify  the  problem  by  using 

xi  -  x.  fi   0  (4) 

where  i  and  j  have  a  common  border. 

This  formulation  can  be  written  in  matrix  form 

Ax^O  (5) 

where  the  elements  of  A_  ,  a. .,  are  either  0,  -1,  or  1.  We  can  now 

enumerate  all  values  of  the  left  side  of  (5) .  The  column  with  the 

_  2 
least  number  of  zeroes  is  now  selected  from  A^^,  a. ,  and  the  value 

of  the  corresponding  x.  is  set  equal  to  one.  The  column  a.  is  then 
added  column  by  column  to  the  matrix  K*^  to  form  a  new  matrix.  The 
process  continues  in  the  same  way  by  searching  the  new  matrix  for  the 
column  with  the  least  number  of  zeros.  In  each  iteration  the  correspond- 
ing x.  is  increased  by  one.  The  column  values  represent  the  value  of 
the  left  side  of  (5) .  A  solution  occurs  when  a  column  is  achieved  that 
has  no  zeros. 


2 
The  speed  with  which  the  algorithm  reaches  a  solution  is  increased 

by  modifying  the  method  of  choosing  a . .  This  modification  is  explained 

in  Section  II.  A.  4.  3 
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The  algorithm  presented  is  based  upon  the  four-color  conjecture, 
but  is  not  restricted  to  the  use  of  only  four  colors,  if  more  colors 
are  desired. 

II.   THE  COLORING  PROCESS 

A.   THE  ALGORITHM 

1.  Let  x.  =  0  for  i  =  1/...,N. 

1 

2.  For  every  2  regions  of  the  map  with  a  common  border  it  is 
required  that  x.  -  x.  ?  0,   i, j ,  =  1, . . . ,N;  i  ^  j. 

3.  Form  the  M  x  N  matrix  (Avbd,)  defined  in  step  2,  where  M  is 
the  number  of  equations  required  to  describe  the  border  constraints 
and  N  is  the  number  of  areas  to  be  colored.  The  elements  of  ^jy^'^-;  > 
are  either  -1,  0,  or  1.  Note: 

a.  For  each  row,  a . ,  i  =  1,. •.,M,  there  are  only  two  non  zero 
elements,  -1,  and  1,  corresponding  to  the  constraints  in  step  2. 

b.  For  column  a.,  j  =  1,...,N,  has  at  least  one  non  zero 
element  and  fewer  than  M  non  zero  elements. 

4.  Let  a.,  j  =  1,...,N,  be  the  number  of  zero  elements  of  each 

a . .  Let  S . ,  j  =  1, . . . ,N,  be  the  number  of  zero  elements  of  each  a. 
%       ,3  3 

which  appear  in  those  positions  where  a . .  =  0 .  Note:  3 .  <  a . , 

j  =  1,...,N.  Find  the  column  of  A^   with  the  smallest  3..  Let  j** 

be  the  smallest  j  such  that  6**  =  min  {£.}.  Let  j*  =  j**. 

:     j    D 

5.  Let  x1  .+  =  x.j.  +  1.  Check  all  the  constraints  in  which  x'  . . 

3*  3*  3* 

appears  to  insure  that  the  constraint  equations  (4)  are  satisfied.  If 

the  constraint  equations  (4)  are  not  satisfied,  continue  increasing 

x.j.  in  increments  of  one  until  either  (4)  is  satisfied  or  x'  ._,,  =  3. 
3  3 

If  x1  .j.  =  3  and  (4)  is  not  satisfied,  set  x'  .#   =  0  and  let  j*  be  the 
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smallest  j  >  j**  such  that  3-*  -  min  (3-)  and  repeat  step  5.  If  no 
such  j  *  exists ,  then  find  the  smallest  a .  and  proceed  as  with  the 
3.'s.  If  no  a .  exists  such  that  x.*  and  x.  satisfy  (4)  for  all  i  7*  j*, 
go  to  step  9. 

6 .  Determine  the  elements  of  column  a . *  by  substituting  the  presently 
defined  values  of  x. ,  i  =  1, — ,N,  into  the  constraint  equations  (4). 
Add  the  column  a.*,  j  =  1,...,N,  to  each  column  of  A^^.  Let  A'^^  = 


A^^  +  [a.  #,...,  a.  *] 


Dx     3' 

7.  Let  a . ,  j  =  1,...,N,  be  the  number  of  zero  elements  of  each  a1 .. 

3  3 

Let  3 . /  j  =  1/ . . . ,N,  be  the  number  of  zero  elements  of  each  a' .  which 

appear  in  those  positions  where  a. .  =  0.  Let  j**  be  the  smallest  j  such 

that  3.  =  min  {3-}  .  Let  j*  =  j**. 
3         j    3 

8 .  Continue  steps  5  through  8  until  a .  =  0 ,  in  which  case  a  solution 
has  been  reached.  Note:  If  x.  =  3  at  any  point  during  steps  5  through 
8,  disregard  the  corresponding  column  a.  in  any  future  iteration,  since 
the  corresponding  x.  is  at  its  maximum  possible  value. 

9.  Let  the  column  a.*  be  the  one  used  to  generate  the  latest  A'  XT 

3*  3  MxN 

matrix.  For  each  zero  element  of  a.*  there  are  2  regions,  g  and  h,  with 
a  common  border  where  x  -  x.  =  0,  i.e. ,  a  violation  of  (4) .  Let  F  be 
the  set  of  all  such  regions  g  and  h.  Denote  by  the  subscript  f  the 
element  of  F. 

10.  For  each  f, increment  the  xf 's  such  that  xf  =  xf  +  1.  Check  all 
of  the  constraints  in  wliich  x^  appears  to  see  if  the  constraint  equations 
(4)  are  satisfied.  If  for  any  f,   an  element  of  F,  the  constraint  equations 
(4)  are  satisfied  set  f  =  j*  and  go  to  step  6.  If  (4)  is  not  satisfied 
for  any  i,  not  an  element  of  F,  where  regions  i  and  f  have  a  common 
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Initial 
AMxN 


^Wn   [ai(D  ^,-ai(l)  ] 


Symbols: 


in  direction  of 
~"^  decreasing  S  .  ,a . 

—  -  >  back  tracking 


AMxN  +  taj(2)  /,,,/aj(2)  ] 

^T 


B3(f 


S^U'+l) 


AMxN+  ^(f)""'^^] 


^ 


^MxN  +  [ajU'+l)  "-'a1(Z'+l)]  3j 


JU) 


'  No  6. 
possible  use 


\       - 

^MxN  +  [aj(£)  '••*,aj(£)] 


(£'+2) 


a 


^MxN  +  [aj  (£+1) 


"1 


3te+i\    i 


\  I 

\ 


^Wn  +  [aj(z'+2)  "'"ai(l'+Z)  ]  ^MxN  +  [aj(&+l)  "-"a-i(i+l)     ] 


3j(0) 


^MxN  +  [a-i(Q)  /-,^ai(Q)3 

Final  solution 
a.  (0)  =  0 


Not  possible  return  to  1**" 
iteration  | 


3J(*+2) 


^MxN  +  [ajU+2)  'm"'aiU+2)l 
Not  possible  return  to 
ft*1  iteration 


Possible  Path  to  a  solution  to  the  map-coloring  problem 

Figure  1 
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border,  note  the  iteration  that  defined  the  present  value  of  x. . 
Repeat  step  10  until  xf  =  3.  Reset  xf  =  0.  Repeat  step  10  for  all 
F.3 

11.  Order  the  iterations  defined  in  step  10  in  decreasing 
iteration  number.  Let  these  iterations  be  called  possible  iterations. 

12.  Reconstruct  the  most  recent  possible  iteration  and  examine 
6 .  =  3  >ilt  or  a .  =  a  .^,  j  =  j*  +  1, . . .  ,N,  depending  on  which  was  used 
to  determine  x.^  in  the  iteration  being  examined.  Choose  a  new  j* 

in  the  same  manner  as  the  j*  was  chosen  in  step  5.  When  a  j*  is  found 
such  that  x.*  and  x.  satisfy  (4)  for  all  i  ^  j*,  go  to  step  6.  If  no 
such  j*  exists  satisfying  (4)  ,  then  the  iteration  being  examined  is  not 
possible.  Repeat  step  12. 

This  completes  the  algorithm.  If  no  possible  solution  exists,  the 
algorithm  will  continue  "backing  up"  the  generated  matrix  chain  until 
the  original  matrix  is  reached.  Figure  1  shows  a  possible  path  to  the 
solution  of  the  map-coloring  problem. 

B.  EXAMPLE 


2    p^l3 

1 

4     J — H 


Fig.  2.  Example  Problem 


The  speed  with  which  the  algorithm  reaches  a  solution  is  increased 
by  starting  the  backtracking  procedure  when  no  variable  associated  with 
the  minimum  6 .  or  a .  can  be  set  at  any  integer  value  between  1  and  3 . 
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where 


xl  "  x2  *  °  x2  -  x3  ^  0 

Xl  "  X3  *   °  X2  "  X4  ^  ° 

xl  "  X4  ^  °  x3  "  x4  *   ° 


set  x.  =  0  for  all  i 

1 

we  write 


xl   x2   X3   X4 
1-10    0 

10-10 

10    0-1 

0    1-10 

0    10-1 

0    0    1-1 


a .    3 

:  


3.   3    3    3    3 
J  


column  used 


then  x,  =  1 

the  column  to  add  (a.^)  is  defined  by: 


xx  -  x2  =  1 

x2  -  x3  =  0 

xx  -  x3  =  1 

x2  -  x4  =  0 

xx  -  x4  =  1 

x3  -  x4  =  0 

where  a.^  =  [111000] 


then  adding  [a .  ^ , . . .  ,a  .*]  +  A^  we  write 
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xl   x2   X3   X4 


0 


0 


0    1-10 
0    10-1 


0    0    1-1 


3     3 


Bi  3 


column  used 


then  x„  =  2 


the  column  to  add  (a.^)  is  defined  by: 


xl  ~  x2    1        x2   x3  "  2 

xx  -  x4  =  1  x3  -  x4  =  0 

where  a.*T  =  [-111220] 
then  adding  [a.*,. . . ,a.  J  +  Aj^  we  write 


xl   xl   x3   X4 


0   -2   -1   -1 


0 


0 

0 

1 

-1 

at 
J 

2 

1 

1 

1 

1 

1 

0 

0 

column  used 
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then  x_  =  3 

the  column  to  add  (a.*)  is  defined  by: 

J 


where  a 


xx  -  x2  =  -1 

xl  "  X3  =  "2 
x±  -   x4  =  1 

=  [-1  -2  1-12  3] 


X2 

"X3 

=  -1 

X2 

"X4 

-  2 

X3 

"X4 

=  3 

then  adding  [a.*, . . .  ,a.J  +  A^^   we  write 


x. 


x. 


x. 


a 


0 

-2 

-1 

■1 

-2 

-3 

2 

1 

1 

•1 

0 

-2 

2 

3 

2 

3 

3 

4 

1 

1 

0 

there  is  no  need  to  calculate  column  a'  because  there  are  no  zeros 
in  column  a'  and  a  solution  has  been  reached.  The  solution  is: 


Area  number 

Color  number 

1 

1 

2 

2 

3 

3 

4 

0 
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III.   CONCLUSION 

The  algorithm  described  will  solve  the  map-ooloring  problem. 
Although  only  four  colors  were  used  in  the  example  in  Section  II-B, 
the  algorithm  can  be  adapted  to  solve  a  problem  for  any  feasible 
number  of  colors.  Several  representative  map-coloring  problems  were 
solved  using  the  enclosed  computer  program.  Appendix  A  shows  a 
solution  to  four- coloring  the  continental  United  States.  In  the 
representative  problems  solved,  both  four  and  five  colors  were  used. 

If  the  solution  is  not  feasible,  e.g.,  the  use  of  four  colors 
when  five  colors  are  required,  as  in  the  split- state  problem,  the 
algorithm  may  consume  several  hours  of  computer  time  before  showing 
that  no  solution  exists.  The  computer  time  is  consumed  in  back- 
tracking the  solution  path  and  in  investigating  every  possible 
solution  path  for  both  the  a ' s  and  the  B ' s . 
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APPENDIX  A 

SAMPLE  PROBLEM  SOLUTION 

The  solution  for  four -coloring  the  continental  United  States  of 
America,  including  boundaries  and  major  water  bodies,  is  as  follows: 

STATE  COLOR 


Maine 

0 

New  Hampshire 
Vermont 

3 
0 

Mas  sachusetts 

2 

Rhode  Island 

3 

Connecticut 

0 

New  York 

3 

New  Jersey 
Pennsylvani  a 
Delaware 

0 
2 
3 

Maryland 
Virginia 
North  Carolina 

0 
2 
3 

South  Carolina 

0 

Georgia 
Florida 

2 
0 

Alabama 

3 

Tennessee 

1 

West  Virginia 
Ohio 

1 
0 

Indiana 

2 

Kentucky 
Illinois 

3 

0 

Wisconsin 

2 

Michigan 
Minnesota 

3 
0 

Iowa 

3 

Missouri 

2 

Arkansas 

3 

Mississippi 
Louisiana 

2 

0 

Texas 

2 

Oklahoma 

0 

Kansas 

3 

Nebraska 

0 

South  Dakota 

1 

North  Dakota 

3 

Montana 

0 

Wyoming 

2 

19 


Colorado 

1 

New  Mexico 

3 

Arizona 

1 

Utah 

3 

Idaho 

1 

Washington 

0 

Oregon 

2 

Nevada 

0 

California 

3 

BORDER  AREAS 

Canada 

2 

Mexico 

0 

Water  -  including: 

1 

Atlantic  Ocean 

Pacific  Ocean 

Great  Lakes 

Great  Salt  Lake 

The  above  problem  was  solved  in  97.2  seconds  on  the  IBM  360/67  computer 
and  used  81  iterations. 
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APPENDIX  B 


PROGRAM  FLOW  CHART 


Assign  initial 
values  to  vari- 
ables 


Read  in  parameter 
and  border 
constraints 


Print  out  answer  f 


Stop 


NZ 


IB 


Yes 


Yes 


Calculate  the  # 
of  common  borders 
for  each  area 


Calculate  aj ,  3  j 
j  =  1,  .^.  ,N  and 
A^  +[aj,...,aj] 


Count  iteration 
number 


I 


Calculate  a , ,x, 
3     J 


Store  a^,x^,j* 
check  X-jft 


Zero  the  counter 

for  a.  &  3- , 
j  =  1,.J.,N  J 


I 


^Wr  +[a-j'---'ai] 


Count  a . ,  3 .  for 
each  column 


I 


Mark  the  column 
with  the  smallest 

3  •  (Denote  as  j*) 


Return 


21 


Calculate  a 


Define  a  dummy 
number  to  indicate 
cannot  use  chosen 


I 


Return 


Return 


Store  aj#,Xj#/j* 


I 


Return 


Yes 


Return  one 
iteration 


22 


Determine  which 
variables  still 
=  0  and  have  com- 
mon  borders 


For  above  vari- 
ables set 


Calculate  which 
iteration  defined 

value  of  x. 

1 


No 


Return 


Define  a  new_ 
column  a,  =  a 


No 


Have 
'all  above 
variables  been 
examined  for  all 
possible  colors? 
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Order  the  iter- 
ations in  de- 
creasing iteration 

# 


Return  to  most 
recent  possible 
iteration 


Assign  to  each 
variable  the  value 
used  prior  to  most 
recent  possible 
iteration 


reconstruct  chosen 
iteration 


Use  next  column 
with  same  number 
of  zeros 


If  no  column  exist 
with  same  number 
of  zeros,  return 
to  next  possible 
iteration 
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& 


No 


Calculate  a.+.x.* 


Store  a.*,x.  ,j* 
-,*  -,*'J 


Return 
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